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ABSTRACT 

We present the ionised gas kinematics of the SBO galaxy NGC 4435 from spectra obtained 
with the Space Telescope Imaging Spectrograph. This galaxy has been selected on the basis 
of its ground-based spectroscopy, for displaying a position velocity diagram consistent with 
the presence of a circumnuclear Keplerian disc rotating around a supermassive black hole 
(SMBH). We obtained the Ha and [N ll] A6583 kinematics kinematics in the galaxy nucleus 
along the major axis and two parallel offset positions. We built a dynamical model of the 
gaseous disc taking into account the whole bidimensional velocity field and the instrumental 
set-up. For the mass of the central SMBH we found an upper limit of 7.5 • 10^ Mq at Scr level. 
This indicates that the mass of SMBH of NGC 4435 is lower than the one expected from the 
M,— (Tc (5 ■ 10^ M0) and near-infrared M,—Lbuige (4 • 10^ M0) relationships. 
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1 INTRODUCTION 

Over the last decade, kinematical studies have proved the presence 
of a supermassive black hole (hereafter SMBH) in the centre of 
about 40 galaxies of different morphological types. For a variety 
of reasons, SMBH are suspected to be present in the centres of all 
galaxies (see Ferrarese & Ford 2005 for a review). 

The census of SMBHs is now large enough to probe the links 
between mass of SMBHs and the global properties of the host 
galaxies. The mass of SMBHs correlates with several properties 
of their hosting spheroid, such as the luminosity (Kormendy & 
Richstone 1995; Magorrian et al. 1998; Marconi & Hunt 2003), the 
central stellar velocity dispersion (Ferrarese & Merritt 2000; Geb- 
hardt et al. 2000), the degree of light concentration (Graham et al. 
2001), and the mass (Haring & Rix 2004). On the other hand, the 
SMBH masses do not correlate with the main properties of discs 
(Kormendy 2001), suggesting that formation of SMBHs is linked 
only to the formation of the spheroidal component of galaxies. The 
recent finding of a new correlation between the central stellar veloc- 
ity dispersion and the rotational circular velocity (Ferrarese 2002; 
Baes et al. 2003; Pizzella et al. 2005) indicates that the mass of the 
SMBH could be also related to the mass of the dark matter halo. 

The mass of SMBHs in elliptical and disc galaxies seems to 
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follow the same scaling relations. However, accurate measurements 
of SMBH masses are available for only 1 1 disc galaxies (Ferrarese 
& Ford 2005) and the addition of new determinations for SO and 
spiral galaxies are highly desirable. Reliable mass estimates of 
SMBHs in disc galaxies have been derived from observations of 
stellar proper motions in the Milky Way (Ghez et al. 2003) and 
spatially resolved kinematics of water vapor masers (Miyoshi et al. 
1995), ionised gas (see Barth 2004 for a review) and stars (see Ko- 
rmendy 2004 for a review) in the other galaxies. It is worth nothing 
that stellar proper motions can only be measured in our galaxy, wa- 
ter masers are not common and both the other two techniques offer 
merits and limitations. Stellar dynamical models are powerful as 
they give information not only on the mass of the SMBH but also 
on the orbital structure of the galaxy. However, both observational 
and computational requirements are expensive. The study of the 
ionised-gas kinematics represents a much easier way to trace the 
gravitational potential of galactic nuclei, since the orbital structure 
of the gas can be assumed to have a simple form, namely the one 
corresponding to a nearly Keplerian velocity field. However, the gas 
is more susceptible to non-gravitational forces and is often found 
in non-equilibrium configurations. Therefore the regularity of the 
ionised gas kinematics has to be verified observationally (Sarzi et 
al. 2001; Ho et al. 2002). The analysis of position-velocity dia- 
grams (hereafter PVD) of the emission-line spectra observed with 
ground-based spectroscopy allows the identification of those galax- 
ies which are good candidates for hosting a circumnuclear Kep- 
lerian disc (hereafter CNKD) rotating around a central mass con- 
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centration (Bertola et al. 1998). Their PVDs are characterized by a 
sharp increase of the velocity gradient toward small radii and the in- 
tensity distribution along the line shows two symmetric peaks with 
respect to the centre (Rubin et al. 1997; Sofue et al. 1998; Funes 
et al. 2002). These objects are good candidates for a spectroscopic 
follow-up with the Hubble Space Telescope (HST). 

In this paper we present and model the ionised-gas kinemat- 
ics of the SBO galaxy NGC 4435 which we measured in spectra 
obtained with the Space Telescope Imaging Spectrograph (STIS). 
This is one of the galaxies we selected on the basis of ground-based 
spectroscopy, for displaying a PVD consistent with the presence of 
a CNKD rotating around a SMBH (Bertola et al. 1998). The paper 
is organized as follows. The criteria of galaxy selection are pre- 
sented in Section|2| STIS observations are described and analyzed 
in Section|3] The resulting ionised-gas kinematics and the morphol- 
ogy of the dust pattern are discussed in Section|4] An upper limit 
for the SMBH mass of NGC 4435 is derived in Section|51 Finally, 
results are discussed in Section^ 



2 SAMPLE SELECTION 

Our sample is constituted by three disc galaxies, namely 
NGC 2179, NGC 4343 and NGC 4435. We considered them as in- 
teresting targets for STIS because ground-based spectroscopic ob- 
servations already allowed the determination of an upper limit for 
their SMBH mass, M, (Bertola et al. 1998). Moreover, they are 
characterized by value of ac~ 150 km s^^. This ctc value is lower 
than those of most of the galaxies so far studied with ionised-gas 
dynamics and higher than the few galaxies studied by means of wa- 
ter masers (Ferrarese & Ford 2005). For this reason, SMBH mass 
determinations in the proposed ctc range would allow a better com- 
parison between data obtained either with gas or stellar dynamics. 
Finally, our sample galaxies belong to morphological types which 
are underrepresented in the sample of galaxies studied so far. 

In this paper we present only the results of NGC 4435, since 
it is the only galaxy of our sample with smooth and circularly sym- 
metric dust lanes as well as it is the only sample galaxy with a 
regular and symmetric velocity field of the ionised gas. This makes 
NGC 4435 an excellent candidate for the dynamical analysis. We 
defer discussion of both NGC 2179 and NGC 4343, as well as the 
estimate of the upper limit of their SMBH mass to a forthcoming 
paper (Corsini et al., in preparation). 

NGC 4435 is a large (2.8 arcmin x 2.0 arcmin, [de Vau- 
couleurs et al. 1991, hereafter RC3]) and bright (Bt = 11.74 
[RC3]) early-type barred galaxy with intermediate inclination (i = 
45°, from RC3 following Guthrie 1992). It is classified SBO°(s) 
and its total absolute mag nitude is M% = -19.41 (RC3), adopting 
a distance of 16 Mpc (Graham et al. 1999). 



3 STIS OBSERVATIONS AND DATA REDUCTION 

The long-slit spectroscopic observations of NGC 4435 were carried 
out with STIS on March 2003 (Prog. Id. 9068, PI. F Bertola). STIS 
mounted the G750M grating centered at Ha with the 0.2 arcsec x 
52 arcsec slit. The detector was the SITe CCD with 1024 x 1024 
pixels of 21 X 21 ^m^. No on-chip binning of the detector pix- 
els yielded a wavelength coverage between about 6290 and 6870 
A with a reciprocal dispersion of 0.554 A pixel"^. The instrumen- 
tal resolution was 1.6 A (FWHM) corresponding to ainstr — 30 
km at Hq. The spatial scale was 0.05071 arcsec pixel"^. 



3.1 Acquisition images 

Four HST orbits were allocated for observing the galaxy. At 
the beginning of the first orbit, two images were taken with the 
F28X20LP long-pass filter to acquire the nucleus. The acquisition 
images have a field of view of 5.4 x 5.4 arcsec^ and a pixel scale of 
0.05071 arcsec. The exposure time was 40 s. The long-pass filter is 
centered at 7230 A and has a FWHM^ 2720 A. It roughly covers 
the R band. 

The first image was obtained by adopting for the nucleus the 
galaxy coordinates from the RC3 catalog. The image was boxcar- 
summed over a check box of 5 x 5 pixels to find the position of 
the intensity peak. The flux-weighted centre of the brightest check 
box was assumed to be coincident with the nucleus location. This 
was used to recenter the nucleus and to obtain the second image. 
After determining the nucleus location, a small move was made 
to centre the nucleus in the slit. The acquisition images were bias- 
subtracted, corrected for hot pixels and cosmic rays, and flat-fielded 
using IRAF^ and the STIS reduction pipeline maintained by the 
Space Telescope Science Institute (Brown et al. 2002). Image align- 
ment and combination were performed using standard tasks in the 
STSDAS package. 

We used the resulting image to check the actual position of the 
slit during the spectroscopic observation as discussed in Section l3!3l 
and shown in Figure^. Moreover, this image was analyzed to map 
the dust distribution in the nuclear region of the galaxy. We con- 
structed the unsharp-masked image using an identical procedure to 
Pizzella et al. (2002). We divided the image by itself after convo- 
lution by a circular Gaussian of width a = 6 pixels, correspond- 
ing to 0.30 arcsec. This technique enhances any surface-brightness 
fluctuation and non-circular structure extending over a spatial re- 
gion comparable to the a of the smoothing Gaussian. The resulting 
unsharp-masked image of NGC 4435 is given in FigureQj. 



3.2 Long-slit spectra 

The position angle of the major axis of NGC 4435 is PA = 13° 
(RC3). We took STIS spectra of NGC4435 with the slit centered on 
the galaxy nucleus and located along its major axis, and with the slit 
parallel to the galaxy major axis on each side of its nucleus with an 
offset of 0.25 arcsec. Following the target acquisition and peak up, 
3 spectra were obtained with the slit along the major axis. Then the 
slit was offset westward by 0.25 arcsec and 4 spectra were obtained 
in the first offset position. Finally, the slit was offset eastward by 
0.5 arcsec and 4 spectra were obtained in the second offset posi- 
tion. Observations of internal line lamps were obtained during each 
orbit for wavelength calibration. At each slit position, subsequent 
spectra were shifted along the slit by 5 detector pixels in order to 
remove the bad pixels. The total exposure time for each slit position 
was balanced within the constraints of a predefined HST offsetting 
pattern. The log of the observations with details about the spectra 
obtained for NGC 4435 are given in TableQ 

The spectra were reduced using the standard STIS reduction 
pipeline. The basic reduction steps included overscan subtraction, 
bias subtraction, dark subtraction, and flatfield correction. Subse- 
quent reduction was performed using standard tasks in the STSDAS 
package of IRAF. Different spectra obtained for the same slit posi- 
tion were aligned using IMSHIFT and knowledge of the adopted 

^ IRAF i.s distributed by NOAO, which is operated by AURA Inc., under 
contract with the National Science Foundation 
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Figure 1. a) HST STIS/F28X20LP acquisition image of NGC 4435. Tlie image lias been rotated to tlie STIS instrumental frame. Orientation and scale are 
given. The white cross is the position of the nucleus from STIS acquisition procedure. The rectangles overplotted on the image show the actual locations of the 
slit during the spectroscopic observations. TablefTllists the offsets of slit position 1, 2, and 3 with respect to the location of the nucleus, b) Unsharp-masking 
of the acquisition image of NGC 4435. c) Portions of the bidimensional STIS spectra of NGC 4435 obtained in position 1 (central panel), 2 (right panel) 
and 3 (left panel). The spectral region centered on the Ho emission line is shown after wavelength calibration, flux calibration and geometrical rectification. 
The spatial axis is horizontal and ranges between —3.2 and +3.2 arcsec, while the wavelength axis is vertical and ranges from 6547 to 6619 A. Individual 
emission lines are identified using a key, such that [N II] (6548), Ha, and [N II] (6583) correspond to [N II] A6548, Ho, and [N II] A6583, respectively, d) Ha 
kinematics from the spectra of NGC 4435 obtained in position 1 (central panels), 2 (right panels) and 3 (left panels). For each slit position the line-of-sight 
velocity curve (top panel), the radial profile of the line-of-sight velocity dispersion (uncorrected for instrumental velocity dispersion, middle panel), and the 
radial profile of line flux in arbitrary units (bottom panel) are given. EiTorbars are not plotted for clarity, but the typical en'ors are shown on the right side of 
each plot, e) Same as in d), but for the [N II] A6583 emission line. 
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Figure 2. Similar to Figure^ but for the archival spectrum of NGC 4435. The acquisition image of NGC 4435 with the slit location, a portion of the two- 
dimensional STIS spectrum after data reduction, and the [N II] A6583 kinematics are shown in panels a), b), and cj, respectively. In a) the image orientation 
and scale is the same as in Figure^. In c) the spatial axis ranges between —3.2 and +3.2 arcsec, and the wavelength axis ranges from 6549 to 6627 A. 



Table 1. Log of the STIS observations of NGC 4435. The number of spec- 
tra, nominal and actual offset, position angle, total exposure time, and ob- 
serving date are listed for each slit position. Slit offsets are given with re- 
spect to the location of the target peak up. 



Pos. 


#Exp. 


Nom. Of. 


Act. Of 


PA 


Exp. T. 


Obs. Date 






(arcsec) 


(ai'csec) 


(°) 


(s) 




1 


3 





+0.05 


12.5 


3012 


9 Mar 2003 


2 


4 


-0.25 


-0.25 


12.5 


3574 


9 Mar 2003 


3 


4 


+0.25 


+0.30 


12.5 


3584 


9 Mar 2003 


4 


3 





-0.08 


89.6 


2673 


26 Apr 1999 



shifts along the slit position. Cosmic ray events and hot pixels 
were removed using the task LACOS.SPEC by van Dokkum (2001). 
Residual bad pixels were corrected by means of a hnear one- 
dimensional interpolation using the data quality files and stacking 
individual spectra with IMCOMBINE. We performed wavelength 
and flux calibration as well as geometrical correction for two- 
dimensional distortion following the standard reduction pipeline 
and applying the X2D task. This task corrected the wavelength scale 
to the heliocentric frame too. The contribution of the sky was deter- 
mined from the edges of the resulting spectra and then subtracted. 
The resulting major-axis and offset spectra of NGC 4435 are plot- 
ted in Figure.^. 

We found in the HST archive 3 other spectra of NGC 4435 
which were taken in 1999 (Prog. Id. 7361, P.I. H.-W. Rix) with the 
same setup we adopted for our observations. These spectra were 
obtained with the slit placed across the galaxy nucleus with a posi- 
tion angle close to the galaxy minor axis (see Table0. We retrieved 
the spectra and reduced them as explained above. The spectrum ob- 
tained close to the minor axis of NGC 4435 is shown in Figure|2j;. 



3.3 Location of the slits 

In our observing strategy, the 0.2— arcsec slit is centered on galaxy 
nucleus and it is aligned along the direction of the columns of the 
acquisition image at the end of target acquisition and peak up. The 
two subsequent offsets were done by applying shifts of —0.25 (i.e. 
westward) and +0.5 arcsec (i.e. eastward) in the direction of the 
rows of the acquisition image. In the archival spectra slit is nomi- 
nally centered on galaxy nucleus. 

We determined the actual location of the slits by comparing 



the light profile of the spectrum with the light profiles extracted 
from the acquisition image. We obtained the light profile of the 
spectrum by collapsing the spectrum along the wavelength direc- 
tion over the spectral range between 6350 and 6800 A. The com- 
parison profiles were extracted from the acquisition image by av- 
eraging 4 adjacent columns. Each strip corresponds to a synthetic 
slit which is 0.2 arcsec wide. Each slit position was determined 
with a minimization of the ratio between the light profile of the 
spectrum and the light profile extracted from the acquisition image. 

We found that the slit centers were misplaced with respect to 
their nominal position. The difference between the nominal and ac- 
tual positions of the slit centre with respect to intensity peak of the 
acquisition image is typically 1 STIS pixel. The locations of the slit 
are overlaid to the acquisition image in Figures.^ and|2^, and the 
details about their positions are listed in TableQ 



3.4 Measurement of the emission Unes 

We derived the kinematics of the ionised-gas component by mea- 
suring the Ha and [N II] A6583 emission lines, which are the 
strongest lines of the observed spectral range. 

For each spectrum we extracted the individual rows out to a 
distance of about 2 arcsec from the slit centre. At larger radii the 
intensity of the emission lines dropped off and we therefore binned 
adjacent spectral rows until a line signal-to-noise ratio S/N ^ 10 
was attained. On each single-row extraction we determined the po- 
sition, FWHM, and flux of the two emission lines by interactively 
fitting one Gaussian to each line plus a straight line to its local 
continuum. The non-linear least-squares minimization was done 
adopting the CURVEFIT routine in IDL^. The centre wavelength of 
the fitting Gaussian was converted into heliocentric velocity in the 
optical convention v — cz. The Gaussian FWHM was converted 
into the velocity dispersion a. The values of heliocentric velocity 
and velocity dispersion include no correction for inclination and 
instrumental velocity dispersion. Errors were calculated by taking 
into account Poisson noise, level of sky and continuum, read-out 
noise and gain of the CCD. The noise associated to the sky level, 
continuum level, and read out was estimated from a spectral region 
free of absorption and emission features after subtracting the fitted 
emission lines. For the archival spectrum of NGC 4435 we fitted 
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only [N II] A6583 line because the Ha line was so weak that little 
kinematic information could be derived from it. 

Heliocentric velocities, velocity dispersions and fluxes mea- 
sured from the Hq and [N II] A6583 lines along the major and offset 
axes of NGC 4435, are plotted in Figures. Ql and^, respectively. 
Heliocentric velocities, velocity dispersions and fluxes measured 
close to the minor axis of NGC 4435 are given in Figure|2j;. 



4 lONISED-GAS KINEMATICS AND DUST 
MORPHOLOGY 

The presence of a fairly-well defined nuclear disc of dust with rel- 
atively smooth and circularly symmetric dust rings and sharply de- 
fined edges is clearly visible in FigureQj. It was first recognized in 
NGC 4435 by Ho et al. (2002). 

We measured the Ha and [N II] A6583 kinematics of 
NGC 4435 out to about 2 arcsec from the centre along the ma- 
jor and offset axes of NGC 4435 (Figures. [T]i and[T];). The Ha 
and [N II] A6583 velocity curves are regular and consistent within 
the errors. However, strong discrepancies are observed along the 
central slit for < r < 0.5 arcsec, and along the western offset 
position for — 1 < r < —0.5 arcsec. We measured consistent val- 
ues of velocity dispersion from the Ha and [N II] A6583 lines, in 
spite of the larger scatter shown by the [N II] A6583 data. 

We measured only the [N II] A6583 kinematics along the slit 
position close to the minor axis of NGC 4435 (Figure|2l;). The mea- 
surements extend out to ~ 1.5 arcsec along the eastern side and to 
~ 0.5 arcsec along the western side, where the dust lanes are more 
prominent. The minor-axis velocity curve and velocity dispersion 
profiles are strongly asymmetric. Our velocities are in agreement 
within the errorbars with those measured by Ho et al. (2002) by 
fitting simultaneously the Ha and the two [N II] emission lines. 

We note that along all slit positions the gas kinematics ob- 
tained from the Ha and [N II] A6583 lines are too different to sim- 
ply combine them in luminosity- weighted mean values. In partic- 
ular, the Ha rotation is characterised by a shallower velocity gra- 
dient than [N II]. Furthermore, most of the Ha flux appear to come 
from circumnuclear regions where low [N II]/Ha ratios suggest star 
formation. On the other hand, the [N II] flux is much more con- 
centrated toward the centre. Since the [N II] emission appears to 
probe better the nuclear regions and is characterised by a simpler 
flux distribution than Ha, we will take into account only the [N II] 
kinematics in our model. 

NGC 4435 has smooth and circularly symmetric dust lanes as 
well as a regular and symmetric velocity field of the ionised gas. On 
the contrary, the other two galaxies of our sample are characterized 
by an irregular kinematics and an irregular dust lane morphology 
(Corsini et al., in preparation). This finding is in agreement with the 
results of Ho et al. (2002), that dust lanes can be used as relatively 
reliable predictor of the regularity of the gas kinematics in the nu- 
clear regions of bulges. This makes NGC 4435 an ideal candidate 
for dynamical modeling. 



5 DYNAMICAL MODEL 
5.1 Velocity field modeling 

5.1.1 Basic steps 

A model of the gas velocity field is generated assuming that the 
ionised-gas component is moving onto circular orbits in an in- 



finitesimally thin disc located in the nucleus of NGC 4435 around 
the SMBH. The model is projected onto the sky plane according to 
the inclination of the gaseous disc. Finally the model is "observed" 
simulating as close as possible the actual setup of STIS spectro- 
scopic observations. The simulated observation depends on width 
and location (namely position angle and offset with respect to the 
centre) of each slit, STIS PSF and charge bleeding between adja- 
cent CCD pixels. The mass of the SMBH is determined by find- 
ing the model parameters which produce the best match to the ob- 
served velocity curve. This modeling technique is similar to that 
adopted by Barth et al. (2001) and Marconi et al. (2003) to analyze 
STIS spectra obtained along parallel positions across the nucleus 
of NGC 3245 and NGC 4041, respectively. 

5.1.2 Model calculation 

Let (r, <j), z) be cylindrical coordinates and consider the gaseous 
disc in the (r, (j>) plane with its centre in the origin. In the case of 
a spherical mass distribution, the gas circular velocity Vc at a given 
radius r is 



Vc{r) = 



GM{r) 



-,1/2 



if) 



2 GM. 

H 



1/2 



(1) 



where M is the total mass enclosed by the circular orbit of radius r, 
(A//L)* is the (constant) mass-to-light ratio of the stellar compo- 
nent (and dark matter halo), and is the circular velocity of radius 
r for a stellar component with {M/L)i,= 1. The radial profile of 
is derived from the observed surface-brightness distribution in 
SectionIO 

The velocity dispersion of the gaseous disc is assumed to be 
isotropic with a Gaussian radial profile 



a{r) — ao + (Tie 



(2) 



Unfortunately, Ha+[N II] imaging at HST resolution is not avail- 
able for NGC 4435. This prevented us to build a map of the surface- 
brightness distribution of the ionised gas, as done for example in 
Barth et al. (2001). We therefore assume that the flux of the gaseous 
disc has an exponential radial profile 

F(r) = Fo + Fie-'-/'-^. (3) 

We now project the velocity field of the gaseous disc on the 
sky plane. Let (x, y, z) be Cartesian coordinates with the origin in 
the centre of the gas disc, the j/-axis aligned along the apparent ma- 
jor axis of the galaxy, and the z-axis along the line of sight directed 
toward the observer The sky plane is confined to the [x, y) plane. 
If the gaseous disc has an inclination angle i (with i = 0° corre- 
sponding to the face-on case), at a given sky point with coordinates 
{x, y) the observed gas velocity v{x, y) is 



■]{x,y) = Vc{x,y)i 



where 



-,1/2 



(4) 

(5) 
(6) 



We assume that the velocity distribution of the gas at position 
(s, y) is a Gaussian with mswv^x, y), dispersion ct{x, y), and area 
F{x,y). 

We now take into account the slit orientation. Let rj, (") be 
Cartesian coordinates with the origin in the STIS focal plane, the 
^— axis aligned with the direction of the slit width, ?7-axis aligned 
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with the direction of the slit length and the (^-axis along the line 
of sight directed toward the observer and crossing the centre of the 
gas disc. The STIS focal plane corresponds to the 77) plane. The 
(a;, y) and rf) coordinate systems are related by the transforma- 
tion 



X = 5 cos 6 — -q sin 6 
y ~ £, sin 6 + r] cos 9 



(7) 

(8) 



where 6 is the angle between the slit direction and the disc major 
axis. 

At position rj) the flux contribution due to gas with a line- 
of-sight velocity v in the range [v — Sv/2, v + 5v/2] is given by 



v+Sv/2 



v~Sv/2 cr(^,r;)v27r 



e.p[V5':)i^.- (9) 



where 5v is the velocity resolution of the model. For a given line- 
of-sight velocity v, F{v) is the "monochromatic" image of the gas 
velocity field observed at A = Ao(l + v/c) for a restframe wave- 
length Aq. 

The LOSVD predicted by the model at position rj) on the 
focal plane is 



(10) 



This takes into account the diffraction of light through the HST and 
STIS aperture. 

For each position rj along the slit the LOSVD predicted by 
the model is given by the contribution of all the points on the focal 
plane inside the slit 



Siv\r,) 



u/2 



(11) 



u/2 



where is the ^ position of the slit centre, w is the slit width, Vd is 
the velocity bin along the wavelength direction in the spectral range 
of interest, and Ma is the anamorphic magnification factor which 
accounts for the different scale in the wavelength and spatial di- 
rection on the focal plane. For STIS in the observed spectral range 
Vd = 25.2 km s~^ pixeP^. The scale in the wavelength direction 
is 0.05477 arcsec pixel" ^ and the scale in the spatial direction is 
0.05071 arcsec pixer \ thus Ma = 0.93 (Bowers & Baum 1998). 
The velocity offset VdMa{£, — Cc) is the shift due to the nonzero 
width of the slit and its projection onto the STIS CCD. The differ- 
ence ^— is in pixel units since Vd is given in km s"^ pixeP^. The 
velocity offset accounts for the fact that the wavelength recorded 
for a photon depends on the position ^ — at which the photon 
enters the slit along the ^ axis (Maciejewki & Binney 2001; Barth 
et al. 2001). This effect is sketched in Figure|3| 

We performed a summation over pixels rather than an integra- 
tion of analytic functions. The model LOSVD and the STIS PSF 
were calculated on a subsampled pixel grid with the bin size of 
5s = 0.01268 X 0.01268 arcsec^ (i.e., a subsampling factor 4x4 
relative to the STIS pixel scale) and on a velocity grid with bin size 
of 5v = 10 km s~^. 5s and 5v correspond to spatial and veloc- 
ity resolution of the model calculation, respectively. In principle, 
smaller values for 5s and 5v could give a more accurate model cal- 
culation, but do not increase the result accuracy. The adopted values 
are the best comprise between good sampling and reasonable com- 
putational time (300 s on a 1-GHz PC). We generated the PSF for 
a monochromatic source at 6600 A using the TINY TIM package in 
IRAF (Krist & Hook 1999). Convolution with the PSF is done using 
the fast Fourier transform algorithm (Press et al. 1992). 
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Figure 3. Schematic representation of the velocity shift caused by the non 
zero width of the slit. The spectra generated along the slit width in pixels 
of coordinates (5i,r)), {^2,v)' UstV)' (54!*?) ^re characterized by a 
velocity offset Av{t;i — ^c), where i = 1,2,3,4 and corresponds to the 
slit centre. The resulting spectrum is S{v\ri). 



The model LOSVD was rebinned on a spatial grid with the 
bin size of 0.05701 x 0.05701 arcsec and on a velocity grid with 
the bin size of Vd ~ 25.2 km s"^ to match the STIS pixel scale 
in the spatial and wavelength direction, respectively. For each slit 
position, the array of model LOVSDs forms a synthetic spectrum 
which is similar to the STIS spectrum. It was convolved with the 
CCD-charge diffusion kernel given by Krist & Hook (1999) in or- 
der to mimic the bleeding of charges between adjacent STIS CCD 
pixels. 

Finally, we analyzed the synthetic spectrum as the STIS spec- 
tra and measured line-of-sight velocity v,nod, velocity dispersion 
amod, and flux Fmod as a function of radius. 



5.2 Stellar component 

To investigate the central mass concentration of NGC 4435 it is 
necessary to determine the contribution of the stellar component to 
the total potential (see Eq.0. 



5.2.1 WFPC2 imaging 

We retrieved Wide Field Planetary Camera 2 (WFPC2) images of 
NGC 4435 from the HST archive. Data for the filter F450W and 
F814W (Prog. Id. 6791, P.I. J. Kenney) were selected to determine 
the stellar light profile minimizing the effects of dust absorption. 
Two images are available for both filters. Total exposure times were 
600 s and 520 s with the F450W and F814W filter, respectively. 
All exposures were taken with the telescope guiding in fine lock, 
which typically gave an rms tracking error of 0.003 arcsec. We 
focused our attention on the Planetary Camera chip (PC) where the 
nucleus of the galaxy was centered for both the filters. This consists 



of 800 X 800 pixels of 0.0455 x 0.0455 arcsec^ each, yielding a 
field of view of about 36 x 36 arcsec^ . 

The images were calibrated using the standard WFPC2 reduc- 
tion pipeline maintained by the Space Telescope Science Institute. 
Reduction steps include bias subtraction, dark current subtraction, 
and flat-fielding are described in detail in Holtzman et al. (1995). 
Subsequent reduction was completed using standard tasks in the 
STSDAS package of IRAF. Bad pixels were corrected by means of 
a linear one-dimensional interpolation using the data quality files 
and the WFIXUP task. Different images of the same filter were 
aligned and combined using IMSHIFT and knowledge of the offset 
shifts. Cosmic ray events were removed using the task CRREJ. The 
cosmic-ray removal and bad pixel correction were checked by in- 
spection of the residual images between the cleaned and combined 
image and each of the original frames. Residual cosmic rays and 
bad pixels in the PC were corrected by manually editing the com- 
bined image with IMEDIT. The sky level (~ 1 count pixeP^) was 
determined from regions free of sources in the Wide Field chips 
and subtracted from the PC frame after appropriate scaling. 

Flux calibration to Vega magnitudes was performed using the 
zero points by Whitmore (1995). To convert to standard B and 
I filters in the Johnson system, we estimated the color correc- 
tion using the SYNPHOT package of IRAF with the SO spectrum 
from Kinney et al. (1996) as template. The color corrections are 
B - F450W = +0.115 and / - F814W = -0.119. 



5.2.2 Correction of dust absorption 

We attempted to correct the data for the effects of dust absorption 
using a method similar to the one described in Cappellari et al. 
(2002). 

For each galaxy pixel we measured the B — I color and we de- 
rived the length a of the semi-major axis of its elliptical isophote. 
The average position angle (6 — 17°) and ellipticity (e = 0.66) 
of the dust features were derived from the analysis of the unsharp- 
masked STIS acquisition images as discussed in Section I531 We 
assumed that the intrinsic galaxy color varies linearly as a function 
of radius. This assumption is justified by Figure|4|which shows the 
intrinsic galaxy color {B — I)q obtained as a straight-line fit to 
the pixel color as a function of a. For each pixel we computed the 
color excess E{B — I) as difference between the measured color 
B — I and the intrinsic galaxy color {B — I)o fitted at that radius. 
This allowed us to obtain a E{V — I) map of the nuclear region 
of NGC 4435. The result is shown in Figure |5| We computed the 
A{I) = 0.558i?(S — /) using the standard Galactic extinction 
curve given by Cardelli et al. (1989) with the assumption that the 
observed color gradient is due to dust rather than stellar population. 
We assumed that dust is distributed in an uniform screen in front of 
the galaxy and we used the A(I) map to correct the /—band image 
for extinction. The extinction-corrected image is shown in Figure 
|5| We applied the A{I) correction only to pixels with \E{B — J)| 
above a given threshold. By comparing the pixel color in the eastern 
and western sides of the galaxy, we defined the threshold as 2 times 
the standard deviation of the observed color with respect to the in- 
trinsic one. We calculated the threshold for 4 < a < 15 arcsec and 
extrapolated it for a ^ 4 arcsec to account for the increasing dust 
absorption at smaller radii (Figure|4}. 

This method corrects the major effects of patchy dust absorp- 
tion. Nevertheless, the dust disc is still visible in the corrected im- 
age, although it is much less optically thick than in the original one. 



Central black hole of NGC 4435 7 



2.01 I ' ' ' I ' ' ' 



1.5 - 




-0.5 - i - 

2 4 6 e 10 IS 14 

a [aresec] 

Figure 4. Calibrated Johnson B — I color for every pixel in the WFPC2 
image of NGC 4435 as a function of its elliptical radius a. Black and red 
dots refer to pixel on the eastern and western side of the galaxy, respectively. 
The thick and two thin continuous lines correspond to titled intrinsic color 
{B — I)o and E{B — I) thresholds, respectively. Thresholds have been 
derived for pixels at 4 < a < 15 arcsec and extrapolated to the innermost 
radii. 



5.2.3 Stellar density profile 

Surface photometry was derived on the /—band extinction- 
corrected image by performing an isophotal analysis with the IRAF 
task ELLIPSE. We derived the isophotal profiles of the galaxy by 
first masking out the remaining dust patches and then fitting el- 
lipses to the isophotes. We allowed the centres of the ellipses to 
vary, to test whether the light distribution in galaxy nucleus was 
still affected by dust obscuration. Since we found some evidence 
of variations in the fitted centre, the ellipse fitting was repeated 
with the ellipse centres fixed to the location found for the outermost 
isophotes. The resulting azimuthally averaged surface-brightness, 
ellipticity, and position angle radial profiles are presented in Figure 
|6| 

The isophotes of the masked image are quite circular with 
e < 0.3 (Figure|6j. For r < 9 pc (2.5 pixels) the ellipticity is poorly 
estimated due to the limited pixel sampling. This allowed us to treat 
the surface-brightness distribution as circularly symmetric, and to 
assume the stellar density distribution as spherically symmetric. 
This approximation is sufficient to estimate the mass-to-light ra- 
tio in the radial range where the ionised-gas kinematics probes the 
galaxy potential. 

We derived the radial profile of the deprojected stellar lumi- 
nosity density T{r) from the radial profile of the observed surface- 
brightness profile E(_R). The intrinsic surface-brightness profile of 
the galaxy S{R) and image PSF were modeled as a sum of Gaus- 
sian components using the Multi Gaussian Expansion (MGE here- 
after) described by Monnet et al. (1992) as done by Sarzi et al. 
(2001). The PSF for the WFPC2/F814W image was generated us- 
ing the TINY TIM package in IRAF (Krist & Hook 1999). The multi- 
Gaussian S{R) was convolved with the multi-Gaussian PSF and 
then compared with E(i?) to obtain optimal scaling coefficients for 
the Gaussian components. The Gaussian width coefficients were 
constrained to be a set of logarithmically spaced values, thus sim- 
plifying the MGE into a general non-negative least-squares prob- 
lem for the corresponding Gaussian amplitudes. For a spherical 
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Figure 5. Dust con-ection of the WFPC2/F814W image of NGC 4435. Left panel: Observed WFPC2/F814W image after calibration to Johnson 7-band. 
Middle panel: Map of the color excess E{B — I). Right panel: Extinction-corrected 7— band image. The field of view is 10.6 X 10.6 arcsec^ and the color 
scale is given in magnitudes. 
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Figure 6. Radial profiles of Johnson /—band surface brightness (top panel), 
ellipticity (middle panel), and position angle (bottom panel) measured on 
the extinction-corrected image of NGC 4435. 
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Figure 7. Deprojection steps for the stellar mass profile of NGC 4435. 
Top panel: Multi-Gaussian fit to observed surface-brightness profile. Middle 
panel: Deprojected stellar luminosity density profile. Bottom panel: Circu- 
lar velocity curve that results if the stellar mass and luminosity density are 
proportional with (M/ L)i,= 1. In each panel the dashed lines correspond 
to the 1(T confidence limits. 



light distribution, the MGE method leads to a straightforward de- 
projection of S{R) into the deprojected stellar luminosity density 
r(r), which can be also expressed as the sum of set of Gaussians. 

For a spherical mass distribution and a radially constant mass- 
to-light ratio, (M/L)*, the stellar mass density p(r) can be ex- 
pressed as p(r) =(M/L)*r(r). It is the sum of spherical mass 
components whose potential can be computed in term of error func- 
tions. The circular velocity u*(r) to be used in Eq.Qis derived as- 
suming {M/L)t= 1. The multi-Gaussian fit to observed surface- 
brightness profile is shown in Figure Q along with the recovered 
luminosity density profile, and the corresponding circular velocity 
curve for (M/L)*= 1. 

We note that since both B and I-band images are affected by 



dust, the B-I maps can only deal a limited description for the dust 
distribution. Unfortunately near-IR images of high spatial resolu- 
tion for NGC 4435 do not exist. The use of such images would 
likely lead to find a larger amount of dust absorption than presently 
estimated. Hence, we are currently underestimating the contribu- 
tion of the stars to the total mass, and overrating the SMBH mass. 
The use of the B-I colour will therefore not affect our conclusions. 
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5.3 Orientation of tlie gaseous disc 

The prediction of the observed gas velocity field for NGC 4435 
depends on the orientation of the nuclear gaseous disc, which can 
be different from that of the main galaxy disc. 

We therefore modeled the gas kinematics by constraining the 
position angle 6 and inclination i of the gaseous disc under the as- 
sumption that the dust lanes are a good tracer of the orientation of 
the gaseous disc (Sarzi et al. 2001). We estimated 8 and i by defin- 
ing ellipses consistent with the morphology of the dust pattern (as 
observed in Figure [ij)) and assuming that dust lanes are circularly 
symmetric. Although the strength of the dust lanes varies spatially, 
it is possible to identify the two most conspicuous ones by visual 
inspection of the unsharp-masked version of the STIS acquisition 
image. To constrain the orientation of the gaseous disc we selected 
two of these features. The inner one has a semi-major axis of ~ 2 
arcsec, which corresponds to the maximum extent of our kinematic 
measurements, and the outer one marks out the edge of dust pattern 
on the western side of the galaxy. 

In each row of the unsharp-masked acquisition image we de- 
termined the position of the two main dust lanes by interactively 
fitting one Gaussian to each absorption feature plus a straight line 
to its local starlight continuum on both side of the nucleus. Once 
determined the position of the dust lanes, we fitted them with two 
ellipses with same centre (Figure |8}. The non-linear least-squares 
minimization was done adopting the CURVEFIT routine in IDL. The 
ellipses have different position angles (8in = 17.0° ± 1.0°, 9out ~ 
13.8° ±0.5°) but same ellipticity(e,„ = eant = 0.66 ±0.03). The 
different position angles of the two ellipses can be interpreted as 
indicative of a slightly warped gaseous disc, but this does not affect 
the final results of the dynamical model as discussed in Section l5!4l 
We assumed for the gaseous disc 6^ = 17° ± 1° and i = 70° ± 2°, 
since the measured kinematics are encircled by the ellipse corre- 
sponding to the inner dust lane. 



5.4 Best-fitting model 

The parameters in our model are the mass M. of the SMBH, the 
mass-to-light ratio (M/L)* of the stellar component (and dark 
matter halo), the inclination i and position angle 9 of the gaseous 
disc, the parameters ao, ai, and of the Gaussian radial profile of 
the intrinsic velocity dispersion of the gas ( see Eq.|2j, and the pa- 
rameters _fo, -Fi, and tf of the exponential radial profile of the gas 
flux ( see Eq.|3}. Given the large number of parameters, it is highly 
desirable to constrain as many as possible of them. We started by 
fixing the orientation of the gaseous disc using the results of the 
analysis of the dust lanes of Section I531 We therefore assumed 
9 = 17° and i — 70° to model the gas velocity field. 

We initially considered a grid of models in which every point 
is determined by AI, and (Af/L)*. For every model in this grid 
we explored several combination for ao, ai, r^, Fo, Fi, and rp 
in order to match the observations. This preliminary analysis re- 
vealed that the predicted flux profile does not depend on the input 
values of M, and (Af/L)*, while the velocity dispersion profile 
depends only on the adopted value of AI,. The flux parameters 
need to be adjusted only if different disc orientations are consid- 
ered. This means that we can adopt the same flux parameters (Fo, 
Fi, and rp) for every point of the grid, but we need to adjust the 
velocity dispersion parameters (ao, cri, ro-) depending on the value 
the black hole mass AI,. Therefore, we can find the optimal values 
for Fo, Fi, and rp at any position in the A4, — (Af/L)* grid, by 
minimizing the Xf = '^{F — Fmod)^ /SF^ where F ± SF and 




Figure 8. The orientation of the gaseous disc of NGC 4435. Diamonds and 
crosses con'esponding to the points of the inner and outer dust lanes are 
overplotted to Figure Absorption features are brighter in this unsharp- 
masked image. The best-fit ellipses of the inner and outer dust lanes and 
their position angles are shown with solid ellipses and solid lines, respec- 
tively. The dashed ellipses and dashed lines enclose the Icr confidence 
ranges on best-fit ellipticity and position angle for both the isophotes. 



F,nod are the observed and the corresponding model flux along the 
major axis (the other slit positions are not reproducible). The best- 
fitting values are Fo = 0.08 and Fi = 0.61 in the adopted arbitrary 
units and rp ~ 10 pc. The velocity dispersion parameters need to 
be optimized for every value of the AI, in the grid, by minimizing 
the Xcr ~ '^{'^ ^ ''"mod)^/5o"^ where a±Sa and a^od are the ob- 
served and the corresponding model velocity dispersion along the 
different slit positions, respectively. The values of the adopted pa- 
rameters for different AI, values are listed in Table|2| Although the 
velocity dispersion parameters depend only mildly on the disc ori- 
entation, we changed them when considering different orientations 
than the one traced by the dust lanes (Section f6.2t . 

We then explored a grid of models with ^ AI, ^ 1.25 ■ 10^ 
Mq and 2.0 {AI/L)-, < 2.35 (M/L)q, where at every point of 
the grid we now use the optimized parameters for the flux and ve- 
locity dispersion obtained before. For each model we calculated the 
= J2i^ ~ ^'mod)^/<5^(f ) where v ± S{v) and v„^ad are the ob- 
served and the corresponding model velocity along the different slit 
positions. The best model has = 955 and a reduced ~ 5.2. 
The best-fitting model requires no AI, and is compared to the ob- 
served [N II] A6583 kinematics in Figure fTTI As we cannot repro- 
duce the observed wiggles and asymmetries, which could be due 
to the patchiness of the surface-brightness distribution of the gas 
(Barth et al. 2001), in a strict x^ sense our best match is not a good 
model. Before deriving confidence limits on AI, and {AI/L)i,, we 
therefore need to allow for the velocity structures that are not repro- 
duced by reseating all x^ values in order for the best to match 
the number of degree of freedom A'^ — AI, where N — 184 is the 
number of the observed points and AI = 2 is the number of the 
parameters in our model. Note that this procedure is equivalent to 
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Figure 9. as a function of M, (/e/f panel) and (M /L)* (ng/if panel). 
The dotted horizontal lines indicate the confidence levels on the best fitting 
values of M,= M© and {M / L)i,= 2.2 (M/L)©. 



allowing for an additional source of error in our data and leads to 
more conservative confidence intervals. Figure|9|show Ict, 2a, and 
3(T confidence levels on M, and {A4/L)t alone, according to the 
Ax'' variations expected for one parameter (i.e. 1, 3, and 9; Press et 
al. 1992). Figure fTol shows similar confidence regions for M, and 
[M/ L)i, jointly, according to the Ax"^ variation expected for two 
parameters (i.e. 2.4, 6.2, 11.8). For the SMBH in NGC 4435 we de- 
rived an upper limit M, ^ 7.5-10'' M© at 3(7 confidence level. The 
mass-to-light ratio is (M/ L)i,— 2.2Q'^'^ at 3(t confidence level. 
The sampling of the grid in the parameters spaces is 2.5 x 10® 
for the SMBH mass, and 0.05 (M/L)© for the mass-to-light ratio. 
The previous results do not change significantly if we use a grid 
with a smaller sampling step. 

As stated before, according to the analysis the model is 
not sophisticated enough to reproduce the data wiggles. However 
to mitigate this problem we can rebin the kinematics in the out- 
ermost regions (beyond \R\ > 0.5 arcsec, where the influence of 
the SMBH is negligible). The errors on every re-binned data were 
computed as the semi difference between the the maximum value 
and the minimum value in the bin. This represent a more reliable 
estimate of the velocity error with respect the error derived from 
the fitting procedure, which take into account the presence of the 
wiggles. We tried several binning on the observed data in order to 
obtain the lower value of the reduced x'^- Using this approach we 
can reach a reduced x^ ~ 1-46, with negligible impact on the pre- 
vious results. The best model still requires no SMBH and the So- 
upper limit is now even tighter (6 - lO" M©), a consequence of the 
augmented relative weight in the of the data close to the cen- 
ter, where the difference between models with different M, is the 
greatest. 

Finally, to account for the uncertainty in the estimate of the 
inclination and position angle of the gaseous disc we calculated the 
best-fitting values of M, and [M/ L)* by building a grid of models 
of the gas velocity fields with 16° sC 6* < 18° and 68° ^ i < 72°. 
We fitted for each disc orientation the velocity dispersion and flux 
radial profiles measured along the major axis. The largest mass for 



Table 2. Parameters of the velocity dispersion radial profile adopted in mod- 
els with (M/L)*= 2.2 (M/L)©, 9 = 17°, i = 70°, Fq = 0.08 and 
Fi = 0.61 in arbitrary units, and rp = 10 pc. 
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Figure 10. grid for models with different values of mass-to-light ratio 
and SMBH mass values. The continuous lines are the confidence levels on 
the best-fitting values of M,= Mq and (M/L)^,= 2.2 (M/L)q . 



SMBH is Af. < 8 - 10® M© at 3o- confidence level to find the opti- 
mized values fov Fq, Fi,rp,aQ, ai, and . The mass-to-light ratio 
ranges between {M/L)^= '^-l^toll and {M/L)^= '2-'25to.lt at 
3a confidence level. These values are consistent within the errors 
with the previously found values. 



6 COMPARISON WITH THE PREDICTIONS OF THE 
NEAR-INFRARED M.-L bulge AND M.-a^ 
RELATIONS 

In order to compare our SMBH mass determination for NGC 4435 
with the predictions of the near-infrared M. — Ltuige relation by 
Marconi & Hunt (2003), we retrieved the 2MASS i/— band images 
from the NASA/IPAC Infrared Science Archive. The galaxy image 
was reduced and flux calibrated with the standard 2MASS extended 
source processor GALWORKS (Jarrett et al. 2000). We chose the 
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Figure 11. Observed [N II] A6583 kinematics (filled circles) along with the best-fitting model (solid line) and model predictions for the SMBH masses of 
NGC 4435 derived in Section l5!4l from the M,—ac relation by Ferrarese & Ford (2005, Af, = 5-10^ Mq, dotted line) and from the near-infrared M,—Lf,y^ig^ 
relation by Marconi et al. (2003, M, = 4 ■ 10^ Mq, dashed line). The observed and modeled velocity curve (top panel), velocity dispersion radial profile 
(middle panel) and flux radial profile (bottom panel) are shown for slit along the major axis (slit #1, upper left panels), close to minor axis (slit #4, upper right 
panel), along eastern offset (slit #3, lower left panel), and western offset (slit #2, lower right panel). 



_ff— band image since it is characterized by lower sky noise with 
respect to the J— and A'— band ones. 

We performed a bidimensional photometric decomposition 
using a Sersic law for the bulge component and an exponential law 
for the disc component and taking into account seeing smearing. 
We adopted the decomposition technique developed by Mendez 
Abreu et al. (2004). The best-fitting parameters are n = 2.39, 
Te = 11.14 arcsec, /i^ — 16.10 mag arcsec^"^, and {b/a)t — 0.64 
for the bulge and h = 14.41 arcsec, /io ~ 18.15 mag arcsec^^, 
and {h/a)d = 0.60 for the disc. The typical error on each parame- 
ter is < 20%. The total bulge luminosity is Lh = 2.8 ■ 10^° Lq, 
which corresponds to a M,— 4-10^ Mq following Marconi & 



Hunt (2003). The presence of the bar in NGC 4435 is enhanced in 
model-subtracted image (Figure fTsl . However, it has to be noticed 
that the bar component dominates the surface-brightness distribu- 
tion of NGC 4435 at a larger radial scale with respect to the exten- 
sion of STIS kinematics. We conclude that the effects of the bar on 
the observed kinematics are negligible in the innermost 2 arcsec. 

Several authors report different values for the central stellar 
velocity dispersion in NGC 4435. 

Bemardi et al. (2002) and Tonry & Davis (1981) both reported 
a value of ac— 174 ± 16 km s~^, whereas Simien & Prugniel 
(1997) measured a value of ac= 156 ± 7 km s~^. The Hypercat 
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Figure 12. Location of the upper limit value of the SMBH mass of 
NGC 4435 with respect to the M, -Uc relation by Feirarese & Ford (2005, 
left panel) and near-infrared A/, — Lj^jg^ relation by Marconi & Hunt 
(2003, right panel). In the left panel, following Ferrarese & Ford (2005) 
we plot the SMBH masses based on resolved dynamical studies of ionised 
gas {open circles), water masers (open squares), and stars {filled circles). 
The upper limit of the SMBH mass of NGC 4335 (Verdoes Kleijn et al. 
2002) has been included for sake of comparison (see Section^. In the right 
panel we plot the SMBH masses which for which /f— band luminosity of 
the host spheroid is available in Marconi & Hunt (2003). In both panels 
dotted lines represent the Icr scatter in M,. 



database also lists an unpublished value of 165 ± 6 km by 
Prugniel & Simien. 

In order to be conservative, we adopted the value (Jc= 156 ± 7 
km s^^ from Simien & Prugniel (1997). Using the correction pro- 
posed by Jorgensen et al. (1995), we derived the velocity dispersion 
that would have been measured within a circular aperture of 1 /8re, 
where is the bulge effective radius derived from our photomet- 
ric decomposition. We found cri/g — 157 kms~^. According to 
the most recent version of the M,—Oc relation (Ferrarese & Ford 
2005), the expected SMBH mass corresponding to this value of 
o"c is 5-10^ Mq . We note that the other published values of ctc would 
lead to even larger black-hole mass for NGC4435. 

In both cases, the upper limit found in our model is signifi- 
cantly below the prediction of these two scaling relations. In Fig- 
ure ll ll we show the comparison between the observed kinematics, 
the best-fit model, and models adopting the prediction of the near- 
infrared M,—Lbuigc (Marconi & Hunt 2003) and M,—Oc (Fer- 
rarese & Ford 2004) relationships, where the intrinsic velocity dis- 
persion profiles were accordingly adjusted. The bigger discrepan- 
cies between the observed and predicted velocity fields are found 
in the innermost ±0.25 arcsec. In Fieure [T2l we report the posi- 
tion of the upper Umit of the M, of NGC 4435 in the M,-ac and 
near-infrared M,—Lbuige relations. 



6.1 Uncertainties due the signal-to-noise ratio 

One important point to investigate is if the S/N ratio of our spectra 
is sufficiently high to exclude the presence of a SMBH with the 
mass predicted by the scaling relations. If it is not the case, the 
central velocity gradient is washed out by the noise and it is not 
possible to detect the presence of a SMBH. 

To exclude this possibility we measured the noise level in the 
observed spectra as the rms of the counts in spectral regions free of 
emission lines. This noise has been added to the model spectra in 
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Figure 14. Comparison between the observed two-dimensional spectra {left 
panels) and the two-dimesional synthetic spectra for (i) the best-fitting 
model, (ii) a model with Al, = 4-10^ Mq as predicted by the near-infrared 
M.-L 

bulge relation, and (iii) a model with Af« — 5 ■ 10^ Mq as predicted 
by the M,—ctc relation for the different slit positions. In the models we 
added a random noise in order to mimic the actual S/N ratio of the ob- 
served spectra. Each box is 15 A X 2.6 arcsec. 



order to mimic the S/N ratio of the observations. Fieure fT4l shows 
that the signature of a SMBH with M, ^4-10^ M© is clearly 
visible in the modeled major-axis spectrum which looks different 
than the observed one. We conclude that the velocity gradient cor- 
responding to M, 4 • 10^ Mq could be measured in the velocity 
profiles of the observed spectra. 



6.2 Uncertainties due to the disc orientation 

The upper limit we derived for A/, depends on the two assumptions 
we did on the orientation of the gaseous disc, i.e. the morphology 
of dust lanes is a good tracer of the orientation of the gaseous disc, 
and the gaseous disc is not warped. 

We test if the presence of a SMBH with a mass of Af.= 4- 10^ 
Mq (i.e., consistent with the near-infrared M.—Ltuige relation) is 
consistent with the major-axis kinematics assuming for the gaseous 
disc a different geometrical configuration with respect to the one 
we derived from the dust lanes analysis. We explored the space of 
parameters 8 and i by building a grid of models of the gas veloc- 
ity fields with M.= 4 • 10^ Mq and {M / L)^,= 2.2. For every 
value of 6 and i we use different values of Fq, Fi, rp, ao, cr-i, and 
Tct, optimising them according to the prescription given in Section 
15.41 in order to match the observations. This value is constrained 
by the gaseous kinematics measured in outer regions which is not 
influenced by the presence of the central SMBH. The best fit to the 
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Figure 13. Two-dimensional pliotometiic decomposition of the surface-brightness distribution of NGC 4435. Left panel: Background- and star-subtracted 
//—band image of NGC 4435 from the 2MASS archive. Central panel: Galaxy model as derived from the photometric decomposition. The bulge is modeled 
vifith a Sersic law (n = 2.39, = 11'.'14, /i,. = 16.10 mag arcsec^^ and b/a = 0.64) the disc is modeled with and exponential law (h = 104'.'41, 
/io = 18.15 mag arcsec"^ and b/a = 0.60). Right panel: Residuals of the model. Images are 3 arcmin X 3 arcmin. 



observed data is found with i = 74° and 9 = 6.5°. Nevertheless, 
in the central ±0.25 arcsec there is a great discrepancy between the 
model and observed kinematics (Figure fTsl which convinced us to 
reject this model as a reliable solution. 

We performed also a simultaneous fit of the velocity curves 
leaving the parameters M, {M/ L)-^, 9 and i free to vary. The best 
fit parameters are M.^ 9 ■ 10*^ Mq, {M/L)-,= 2.5 ± 2 (M/L)© 
9 = 15° ± 2°, and i = 75° ± 3°. Although for computational 
reasons this minimization is performed with all parameters for the 
intrinsic flux and velocity dispersion profiles fixed at the best values 
derived in Section l?!4l the similarity of the disc orientation found 
in this fit with respect to the one determined in Section l5!3l makes 
us confident about the choice of adopting the dust-lane morphology 
to constrain the orientation of the gaseous disc. 

The previous results still rely on the assumption that the gas 
follows a coplanar distribution. However, the presence of a SMBH 
with a mass Af,= 4 ■ 10^ Mq could still be consistent with ob- 
served kinematics if the gas would warp to a more face-on orienta- 
tion toward the central region. By exploring the 9 and i parameter 
space with models for the gas velocity field only within ±0.25 arc- 
sec from the centre and with M,= 4 ■ 10^ M© and {M/L)^,= 2.2 
(M/L)o, we find that the observed kinematics could be explained 
if i = 28° and 9 = 10° (Figure [T^l. The images, however, do 
not suggest such dramatic variation of the disc orientation, and in 
Figure[8|in particular the dust lanes suggest an highly inclined con- 
figuration down to very small radii (r ~ 0.5 arcsec). 



6.3 Uncertainties due to the slit position 

The upper limit we derived depends also on the assumed position 
of the slits. In Section l33l we determined the actual position of the 
slits by comparing the light profile of the spectrum with the light 
profiles extracted from the acquisition image. In the case that this 
procedure is not reliable, we try to find new slits position (but not 
changing their relative distance) in order to allow the presence of a 
M,= 4 • lO'^ Mq. We used the parameters 9 = 17° and i = 70° 
and (M / L)-t,= 2.2 (M/L)© and we find that the major axis velocity 




K [arcsec] 



Figure 15. Observed [N II] A6583 (filled circles) and fitted (solid line) kine- 
matics in the innermost ±2.5 arcsec along the major axis of NGC 4435 
obtained with M. = 4 • 10^ Mq and (M/L)t= 2.2 (M/L)q but allowing 
the position angle and inclination of the gaseous disc to vary. The central 
velocity gradient is not reproduced by any choice of disc inclination and 
position angle. 



gradient is well reproduced if we shift the slit of 3 S T I S pixels. But 
in that case, the kinematics along the parallel offset slits is not re- 
produced anymore: there are no ways to reproduce simultaneously 
all the slits. Having parallell slits provide a kinematic test for cor- 
rect slit positioning, which confirmed the slit locations inferred in 
Section l3!3l using the acquisition images. 
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Figure 16. Observed [N II] A6583 (filled circles) and fitted (solid and 
dotted lines) kinematics in tiie inneitnost ±2.5 arcsec along the major 
axis of NGC 4435. The solid lines correspond to a model obtained with 
M,= 4 • 10^ M0, {M/L)i,= 2.2 (MIL)q, 9 = 10° and i = 28° in the 
innermost ±2.5 arcsec. The dashed lines correspond to a model obtained 
with M,= Mq, (M/L)i,= 2.2 (M/L)©, 9 = 17° and i = 70°. 
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Figure 17. Observed [N II] A6583 (filled circles) and fitted (solid lines) 
kinematics in the innermost ±2.5 arcsec along the major axis of NGC 4435. 
The model is obtained using a M.=4 ■ lO'' M@i = 70°, 9 = 17°, 
(M/L)i,= 2.2 (M/L)0, and an ad hoc flux distribution to reproduce the 
central gradient. The predicted flux radial profile is completely inconsistent 
with the observed one. 



6.4 Effects of the flux distribution on the measurements 

Without narrow-band images we had to resort an analytical de- 
scription for the intrinsic surface brightness of the ionised gas 
(see Section l5.1.2t . Our choice of an exponential radial profile for 
the emission-line fluxes seems well-justified a posteriori given the 
good match to observed flux profiles along all slit positions, partic- 
ularly towards the centre (see Figure fTTl . However, as in the case of 
NGC 3245 (Barth et al. 2001), it is likely that unaccounted small- 
scale flux fluctuations are responsible for the poor match in many 
part of the velocity curve (see Section l5!4l 

As far as the SMBH mass measurement is concerned, it is im- 
portant to bear in mind our ignorance of the details characterising 
the gas surface brightness towards the very centre of the NGC 4435. 
In particular, the presence of nuclear dust may represent a limitation 
for our analysis. If the very central regions, where the gas clouds are 
more directly affected by the gravitational pull of the SMBH, are 
affected by significant dust absorption, it possible that gas further 
away from the centre and moving at slower speed could contribute 
to the observed central velocity gradient and flux peak. If this is the 
case we would be underestimating the SMBH mass. 

We investigated this possibility by using a different model for 
the intrinsic gas surface-brightness, including an additional central 
component with negative amplitude, to mimic a central flux deple- 
tion due to dust absorption. We have started by trying to accomo- 
date a 4 • lO'^ M0 SMBH. As shown by Figure fTTl to match the 
observed gas rotation significant dust absorption is required out to 
about 1 arcsec, but in this case the predicted flux distribution is in- 
consistent with the observed profile. The extent of the nuclear dust 
is a direct consequence of the fact that the observed rotation curve 
can be fitted without a SMBH (Figure fTTl , so that dust have to 
screen the emission from regions where the SMBH can contribute 
significantly to the observed rotation curve. For a (M/L)*= 2.2, 



including a 4 ■ 10^ M© SMBH still increases the circular velocity 
(Eq. by 10 km at a radius of 55 pc, corresponding to 0.7 
arcsec. 

In this framework, to explain the observed flux profile we 
would need to resort to fine-tuned geometries for the gas distri- 
bution in which intervening off-plane material would contribute 
to the observed central peak in flux distribution. However, the ex- 
ceptional regularity of the dust distribution in our images and the 
well-established connection between dust and gas (e.g., Pogge et 
al. 2000) strongly suggest that the gas resides on a simple a simple 
disc. 

Furthermore, the required amount of central flux depletion 
needed to accomodate a 4 ■ 10*^ Mq SMBH would translate into 
an additional central reddening of E{B — I) — 1.3 mag. which is 
should be easily detected in our map for the E{B — I) color ex- 
cess (Figure 5, in which the central value is 0.6 mag.). This is not 
the case, however. In this calculation we had taken into account i) 
the different value of the extinction Ax at different wavelenght (the 
J— band and the [N II] A6583 region , using the standard Galactic 
extinction as done in Cardelli et al. 1989); ii) that only the light 
from stars behind the dusty disc is absorbed. 

We therefore consider it unlikely that nuclear dust could hide 
a SMBH consistent with the expectations of the M,—ac relation. 
On the other hand, it is still possible that nuclear dust could impact 
our inferred upper limit on the SMBH mass while still matching 
the flux profile. In this case however, the central dust cannot extend 
beyond the width of our slit. The impact on Al, will therefore be 
rather limited; we find that our upper limit would increase by only 
few percent, to 8 • l(f Mq . Gaussian profile instead of an exponen- 
tial profile for the flux distribution, as done by Sarzi et al. (2001, 
2002), would have a similar small effect. 
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6.5 Correction for asymmetric drift 

The velocity dispersion of tlie gas peaks at 100 km s^^ in the cen- 
tre. This intrinsic line width is far higher than that expected ei- 
ther from rotational and instrumental broadening or from thermal 
motion. This finding is in agreement with earlier results on bulges 
of other disc galaxies based on both ground-based (Fillmore et al. 
1986; Bertolaetal. 1995; Cinzano etal. 1999; Pignatelh et al. 2001) 
and STIS spectroscopy (Barth et al. 2001) and suggests that ran- 
dom motions may be crucial for the dynamical support of the gas. 
If ionised-gas disc is fragmented into coUisionless clouds, then dy- 
namical pressure supports it against gravity and the mean rotational 
velocity is smaller than the circular velocity Vc given in Equa- 
tionQ As a consequence, the enclosed mass is underestimated by 
dynamical models which do not account for the asymmetric drift 
correction. 

In the following we address how much of an effect this would 
have on the difference between the velocity gradient we measured 
and that we inferred for gas rotating around a SMBH with M,= 
4 • 10^ Mq as predicted by the near-infrared M,—L),uige scaling 
relation. 

We assume that the motions in the gas disc are close to 
isotropic in the radial and vertical direction. Then ~ <Jr and 
(vrVz) = 0, and the asymmetric drift correction can be expressed 
as 
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where cr^ is the azimuthal velocity dispersion and u is the number 
density of gas clouds in the disc (e.g., Binney & Tremaine 1987). 
We assume that the number density of the gas clouds is propor- 
tional to the gas flux whose radial profile has been parametrized 
in Eq.|3| The relation between Or and is given by the epicycle 
approximation 
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The approximation is strictly valid in the limit a <^ Wc. If we as- 
sume that the SMBH mass is given by the near-infrared M,—Lf,uige 
relation, we have a ratio (Jr/v^ < 0.32 everywhere in the gaseous 
disc of NGC 4435 (Figure fTsl and the method can be adopted for at 
least an approximate treatment of the asymmetric drift. After calcu- 
lating the asymmetric drift, we obtained the line-of-sight velocity 
and velocity dispersion as 



and 
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respectively. These values of v and a are used in Eq. |9| for the 
model calculation. For models without an asymmetric drift, the 
above equations reduce to Vc = and a — ar = u^, correspond- 
ing to the case where the intrinsic velocity dispersion represents a 
locally isotropic turbulence within the gaseous disc. 

We applied the asymmetric drift correction for M,= 4 • 10^ 
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Figure 18. The ratio <Jr/vc for the model with M,= 4 ■ 10^ Mq. 
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Figure 19. Upper panel: The observed [N II] A6583 (filled circles) along 
the the major axis of NGC 4435 is compared to the one obtained with Mm = 
4 ■ 10'^ Mq, {M/L)i,= 2.2 (M/L)q, 9 = 17° and i = 70° with (solid 
line) and without (dashed line) applying the asymmetric drift correction. 
Lower panel: Ratio between the rotation curves obtained with and without 
applying the asymmetric drift correction. 



Mq and found that the contribution of intrinsic velocity dispersion 
to the rotation curve is negligible (^ 5%). The comparison between 
the model rotation curves obtained along the galaxy major axis with 
and without the asymmetric drift correction is shown in Figure fT9l 
We can conclude that the contribution of the velocity dispersion is 
not enough to allow the presence of M,— 4- 10^ Mq in the nucleus 
of NGC 4435. 



7 CONCLUSIONS 

We presented long-slit STIS measurements of the ionised-gas kine- 
matics in the nucleus of the SBO NGC 4435. NGC 4435 belongs to 
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a sample of 3 galaxies which were selected as STIS targets on the 
basis of their CNKD kinematics (Bertola et al. 1998). 
For NGC 4435 we found the following: 

(i) It shows a CNKD with a symmetric and regular ionised-gas 
kinematics, which is suitable for dynamical modeling in order to 
estimate the mass of the central SMBH. Moreover, NGC 4435 is 
characterized by the presence of smooth and circularly symmetric 
dust lanes. This result is in agreement with early findings of Ho 
et al. (2002) who found an empirical correlation between the mor- 
phology of the dust lanes and the degree of regularity in the gas 
velocity field. 

Bertola et al. (1998) demonstrated that it is possible to de- 
tect the signature of a CNKD in the emission-line PVDs obtained 
from ground-based spectroscopy of nearby galaxies. Using this re- 
sult, Funes et al. (2002) estimated that the frequency of CNKDs 
in randomly selected emission-line disc galaxies is < 20%. How- 
ever, this criterion has to be combined with that of the presence of 
a regular dust-lane morphology (Ho et al. 2002). On the basis of 
our STIS observations, we conclude that less than 10% of nearby 
disc galaxies, which show narrow emission lines in their ground- 
based spectra, host a CNKD with a velocity field which is useful 
for dynamical modeling at HST resolution. 

(ii) We modeled the ionised-gas velocity field of NGC 4435 by as- 
suming the gas is moving onto circular orbits in an infinitesimally 
thin disc located around the SMBH. We constrained the orientation 
of the gaseous disc (6 = 17° ± 1°, i = 70° ± 2°) by assum- 
ing that the morphology of dust lanes is tracing the gas distribu- 
tion. The comparison between the observed and modeled velocity 
field has been done by taking into account several effects related 
to telescope and instrument optics (nonzero aperture size, apparent 
wavelength shifts for light entering the slit off centre, and anamor- 
phic magnification) as well as detector readout (charge bleeding 
between adjacent CCD pixels). We derived for the SMBH mass of 
NGC 4435 an upper limit M. ^ 7.5 • 10*^ M© at 3cr confidence 
level. 

(Hi) The upper limit we derived for the SMBH mass of NGC 4435 
considerably falls short of the value of M,— 4 ■ 10^ Mq pre- 
dicted by the M.— tJc relation (Ferrarese & Ford 2005), as well 
as of the value of M,— 5 • 10^ Mq predicted by the near-infrared 
M,—Lbuige relation (Marconi & Hunt 2003) as shown in Figure 
1121 This discrepancy is not due to the noise in our spectra (Section 

6.1) , is robust against uncertainties in the disc orientation (Section 

6.2) or slit positioning (Section 6.3), and cannot be ascribed to the 
presence of nuclear dust (Section 6.4) or pressure-supported gas 
(Section 6.5). The measured velocity gradient is consistent with 
a M.= 4 • lO'^ M© only if we assume an ad hoc orientation 
(6^ = 10° ± 1°, i = 28° ± 2°) for the central portion of the gaseous 
disc (|r| ^ 0.25 arcsec). However, the images do not suggest such 
a dramatic variation of the disc orientation in the observed radial 
range. 

(iv) The case of NGC 4435 is very similar to that of the ellipti- 
cal NGC 4335 (Verdoes Kleijn et al. 2002), where the ionised-gas 
kinematics measured with STIS along three parallel positions is 
consistent with that of a CNKD rotating around a SMBH with an 
unusually low mass {M, ^5 10* Mq) for its high velocity disper- 
sion ((Tc= 300 ± 25 km see FigurefTH. Both NGC 4435 and 
NGC 4335 are amenable to gas dynamical modeling. Their nuclear 
discs of ionised gas have a well-constrained orientation as well as 
a symmetric and regular velocity field. There is no evidence that 
gas is affected by non-gravitational motions or that it is not in equi- 



librium, as observed for instance in the E3 galaxy IC 1459 (Cap- 
pellari et al. 2002) and the Sbc NGC 4041 (Marconi et al. 2003). 
For NGC 4335 the observed high gas velocity dispersion in the 
central arcsec potentially invalidates the M, measurement based 
on a thin rotating disc model leading to a larger M, value and to 
a smaller discrepancy with the M, scaling laws. On the contrary, 
in NGC 4435 the more limited width of the lines allowed us to 
consider the possibility of a modest dynamical support for the gas 
clouds. This has little impact on the measured M, and can not ex- 
plain its discrepancy with the values predicted by the M, scaling 
relations. As such, if we suppose the resulting upper limit for the 
SMBH mass of NGC 4435 is also unreliable, we are then forced to 
conclude that all the SMBH masses derived from gaseous kinemat- 
ics following standard assumptions have to be treated with caution, 
even if they agree with the predictions of the scaling relations. 

Thus, NGC 4435 is so far the best candidate of a galaxy with a 
massive bulge component with a lower M, content than predicted 
by the M,—cjc and Mt—Lbuige scaling relations. Nevertheless, 
there is an increasing evidence of galaxies for which dynamical 
models impose an upper limit to the SMBH mass which is either 
lower than or marginally consistent with the one predicted by the 
M.-(Jc (Sarzi et al. 2001; Merritt et al. 2001; Gebhardt et al. 2001; 
Valluri et al. 2005). These objects may represent the population of 
laggard SMBHs with masses below the M,—ac relation as dis- 
cussed by Vittorini et al. (2005). Laggard SMBHs are expected to 
reside in galaxies that spent most of their lifetime in the field, where 
encounters are late and rare so as to cause only slow gas fueling of 
the galactic centre and a limited growth of the SMBH mass. 

In light of the serious implications on the slope and intrin- 
sic scatter of the A/, relations that our results would imply, an in- 
dependent measurement of the SMBH mass NGC 4435 based on 
stellar dynamical modeling is highly desirable. This could be done 
with near-infrared spectroscopic observations with 8m-class tele- 
scopes assisted by adaptive optics systems. Recently, Houghton et 
al. (2005) have shown in the case of the giant elliptical NGC 1399 
that this technique can deliver diffraction limited high signal-to- 
noise spectra suitable for measuring the stellar LOSVD within the 
sphere-of-influence of the SMBH, providing a concrete alternative 
to HST that today is in serious jeopardy. 
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APPENDIX A: CALCULATION OF THE COLOR EXCESS 
E{B - I) IN THE CASE OF A M.= 4 • 10^ M© 

In this section we evaluate the colour excess E'{B ~ I) we would 
expect in the case that a M, = 4-10^ M© is present and the absorp- 
tion of the central dust is able to reproduce the observed kinematics 
(see Section 6.4). 

If the ionised gas is settled in the equatorial plane as well as 
the dust, the intrinsic flux corrected by absorption (f will be 



related to the intrinsic one (F""^") by the relation 
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The coefficient c indicates the efficiency of the absorption, and it 
can be easily derived (see Figure lATl . For c = 1, F'^'"' = 0, for 

^ _ Q pahs _ perns 

Using the standard Galactic extinction (Cardelli et al. 1989) 
we can determine the relation between the extinction A[jv7/] and 
the extinction in the I band, Aj: 

Al = 0.588yl[iv„] (A2) 

The central stellar component (bulge plus disc) is thicker with re- 
spect the ionised gas and the dust discs. Thus, only the portion of 
stars behind the dust layer is absorbed: 
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What we should observe is the extinction A'^ defined by: 

f:^^ = io-°'*-"'^ ■ (A4) 



with (comparing Eqs. A3 and A4): 
1 -I- in-" * '*/ 



(A5) 



Then, according to the standard Galactic extinction (Cardelli et al. 
1989), the expected color excess E'{B — I) is: 



E'{B - /) = A'j/0.5583 



(A6) 



Then, combining Eqs. Al, A2, A5 and A6, it is possible to calculate 
the map of E'{B — /) as a function of c, derived from Figure lATI 

1 I (-1 „^0.588 

E'{B -I) = -4.48 log \ (A7) 

In the central region, where c = 1 the expected color excess is 1.3 
mag, while the observed one is 0.6 mag. (Figure 5). 



This paper has been typeset from a T^f/ ETbX file prepared by the 
author. 
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Figure Al. Intrinsic flux radial profile of the gas in the NGC 4435 model. 
The solid thin line shows the intrinsic emission of our best model, whereas 
the solid thick line shows the absorbed emission needed to accomodate a 
M,= 4 • lO''^ Mq and match the observed velocities (Figure 17). The 
dashed line shows the efficiency of the dust absorption, c, as defined in 
the text. 



